Narrowband frequency tunable light source of continuous quadrature entanglement 
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We report the observation of non-classical quantum correlations of continuous light variables from 
a novel type of source. It is a frequency non-degenerate optical parametric oscillator below threshold, 
where signal and idler fields are separated by 740MHz corresponding to two free spectrum ranges of 
the parametric oscillator cavity. The degree of entanglement observed, - 3.8 dB, is the highest to- 
date for a narrowband tunable source suitable for atomic quantum memory and other applications 
in atomic physics. Finally we use the latter to visualize the Einstein-Podolsky-Rosen paradox. 



(N 
O 

o 



(N 
> 

in . 

o ■ 
in : 
o . 

(N ■ 

o ; 
^ ■ 

Oh! 
+- > ■ 

C ■ 
d • 

s : 
cr 



X 
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I. INTRODUCTION 



Quantum interface between light and atomic systems 
requires narrowband frequency tunable sources of light 
possessing various non-classical features. Examples of 
sources which are in demand include narrowband tunable 
single photon sources and narrowband tunable sources of 
Einstein-Podolsky-Rosen entangled light. Both types of 
sources can be constructed using a non-degenerate op- 
tical parametric oscillator below threshold (ND-OPO). 
However, the vast majority of current experiments with 
entangled light is carried out with single pass downcon- 
version sources with their inherent low spectral bright- 
ness. The first EPR-type experiment using a narrow- 
band light source jj], ^ was also the first demonstration 
of entanglement for continuous variables. Recently an- 
other experiment demonstrating the EPR state for soli- 
tons has been carried out || although in this experiment 
the bandwidth of the correlated fields was very wide. The 
first demonstration of a frequency tunable source of a 
non-classical narrowband light source possessing entan- 
glement between the modes at 852 nm and 917 nm was 
reported in Q . A potentially tunable source of entangled 
beams at the same wavelength but in different spatial 
modes was used in the continuous variables teleporta- 
tion experiment J5j. In the latter experiment a degen- 
erate OPO pumped in two opposite directions was uti- 
lized to generate two squeezed beams which were mixed 
on a beamsplitter to produce two entangled beams with 
the degree of entanglement around -1 dB. In the present 
paper we demonstrate the entanglement of a new type. 
Namely two symmetric longitudinal modes of the OPO 
cavity which are shown to possess the quadrature phase 
entanglement as suggested in the original paper by M. 
Reid where entanglement of the OPO output has been 
discussed for the first time. 



II. ND-OPO THEORY 

In this section we present a model for calculating quan- 
tum fluctuations of the ND-OPO output modes. Our 
treatment is based on || with experimentally relevant 
generalizations: the two modes do not necessarily expe- 
rience the same losses in the resonator and the output 
coupling rates can be different too. The ND-OPO is de- 
scribed as a moderate finesse optical resonator containing 
a x^ 2 -* nonlinear medium. This resonator is assumed to 
be resonant with two field modes a+ and a_ , henceforth 
denoted the signal and idler modes respectively. In terms 
of frequency the signal and idler modes are assumed to 
be an integer number of free spectrum ranges higher and 
lower relative to the degenerate frequency set by the clas- 
sical external pump field a. That is, energy conservation 
requires that the sum of the frequencies of the resonating 
fields must equal the pump field frequency cj + +cj_ = u p , 
so that lu± — lUp/2±vujfsr, where Ufsr is the frequency 
of a free spectrum range and v is an integer. The mode 
operators obey the usual bosonic commutation relation 

\a n (t), a^(t)] = S nm , {n,m} = {+,-} (1) 



A. Dynamics 

The starting point for our treatment of the ND-OPO 
is the Hamiltonian 



H = H 



free 



(2) 



where free, int and res describes the free, noninteract- 
ing OPO cavity fields, the nonlinear interaction between 
the fields and the coupling of the fields to the external 
reservoir of vacuum modes respectively. The first two 
Hamiltonians take the form 



H 



free 



(3) 
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and 



Hint = ihg {ae lx a\ 



at - a*e lx a + a_ ) (4) 
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where g is the coupling strength, a = a exp (—2iujt) is the 
classical pump field and x is the pump field phase. From 
(||), (0) and (^) we can derive the Heisenberg-Langevin 
equations of motion for the cavity fields 

■£a± = -iu±a± + gae^a^ - j±a ± - , 

+V^SS l + V2(7±-k±)6£ 

where the a m fields describe the vacuum fluctuations 
leaking into the resonator through the main coupler and 
the b ln fields enter the ND-OPO via residual leakage. 
k± define the leak rate through the main coupler at fre- 
quencies lu± and 7± define the overall leak rates. Our 
observables of interest are the slowly varying envelopes 
of the field quadratures 

q±(0) = a ± e* u ±*- e ) +4e" i(w±t_e) (6) 

The pump field we parameterize in terms of its threshold 
value a = sctth, where ath — \J 1+1-/9- This results in 
the equations of motion for the subharmonic quadrature 
operators 

f t q±{6) = syiTl^QT ix - 0) - l± q±(0) (7) 

+V^q±(8) + V2( 7 ±-K±y± n W 1 ' 

Here we have defined the slowly varying vacuum field 
quadratures q£(0) = a^e^^-^ + ia^e' 1 ^-^ and 



p™(0) = l™ e i( u± t-e) + (p ± yn e -i(u, ± t-e) ^ The phage dg _ 

pendence of q ln and p m will be left out in the following 
since no phase can be assigned to the vacuum state and 
hence the properties of the 'in' quadrature phase opera- 
tors is independent of the choice of 6. 

From (0) we see that only the quadratures with phases 
4>± = x/2 ± 8 will interact. Hence we find that the ND- 
OPO intracavity fields are described by a continuum of 
pairwise coupled equations of motion of the form 



3S?±(#±) = £ \/1+1-9t (^ t ) ~ 7±g±0 ±) 



Before moving on to solve (^) we make a few observa- 
tions about the quadrature phase operators. Using (|l|) 
we find that these operators obey the following commu- 
tation relations 

[Qn(<j>), qm(4>')} = 2i sin(<f>' - <j))6 nm , {m, n} = {+, -} 

meaning that the combined signal/idler quadratures de- 
scribing cross correlations of the ND-OPO output com- 
mute as 



+ V)] = [q+(<f>+)-q- 
I 



+ ip)) = 4zsin-0 (10) 



Hence the combinations {q + ((j) + ) + g_ ((/>_), q+{<p+ + 
it/2) + q-{4>- + tt/2)} and {q+{<p+) - q~ (</>-), <7+(0+ + 
it/2) — q-((j>- +7r/2)} form conjugate observables. As we 
will see below the fluctuations of the quadrature differ- 
ence q+{<f)+) — q— ((ft—) can become very small in the ND- 
OPO output as compared to a two mode coherent state 
output. This means that +7r/2) + g_(</)_ +7r/2) = 

q+W + ) - ?_(#_), with <f>' ± = x/2 ±(0 + n/2), also 
can become a quiet observable since <j>± is defined for 
any arbitrary value of 9. From the commutation rela- 
tion (10) we then conclude that q + ((f> + ) + g_(0_) and 
q+(<j>+ + tt/2) — q_(0_ + tt/2) must become correspond- 
ingly noisy observables. This is illustrated in the phasor 
diagram of Fig. [I} In the parametric down conversion 
the phase is conserved, meaning that the signal and idler 
fields are created symmetrically around the axis making 
an angle of x/2 with the abscissa in phase space. This 
axis illustrates the pump phase. In the rotating frame 
oscillating at ui p /2 the signal field will rotate clockwise 
and the idler will rotate counterclockwise in such a way 
that their sum always will be parallel to the pump axis. 
Hence we find that the difference of projections of sig- 



nal/idler fields onto the pump axis will tend to cancel 
on average. This difference is described by the cross 
correlation of quadratures </+(</>+) — q_(0_). Also the 
sum of projections onto the axis perpendicular to the 
pump axis will cancel. This observable is described by 
+ tt/2) + q-(cf>-+ tt/2). 



B. Quantum fluctuations 

The solutions of @ are found by Fourier transforming 
and solving the resulting algebraic equations 



Vq = Aq tn +Bf 
where we have defined 



(11) 



V 



_( J+-iti -£^1+1^ \ a= (V2k+~ I) 



-£^7+7- 7- - ^ 
B 



^/2~^Z 







^2(7+ - k+) 

V 2 (7- - k-) 



(12) 
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FIG. 1: The parametric downconversion seen in a phasor 
diagram. In a frame rotating at the degeneracy frequency a+ 
is rotating clockwise and a_ is rotating counterclockwise with 
relative phases such that the two phasors always add up to 
be parallel to the dashed line. 

and 

<-(t\£l) •<*-($)•''-($) < 13 > 

now are the Fourier transformed quadrature phase oper- 
ators. From the solution of ([ll]) we can find the fields 
escaping the ND-OPO through the main coupler via the 
boundary condition on the coupler 

«f »* = Aq- T^ n (14) 

These are given by 

q° ut = [AD- 1 A - 2] q^ n + AD~ x Bp in (15) 

where X is the 2x2 unit matrix. To describe the non- 
classical correlations between the out fields in ([Tot), we 



define the general linear combinations of the signal and 
idler quadratures 

Q=( % + )=Uq, U=( C0S ^, Shi t) (16) 
\Q- J K—smtf) cost/" J v ' 

Hence we find the outgoing linear combinations as 

Q out =U [AV^A-l^Q" 1 +UAV- 1 BU T P in (17) 

where Q m = Uq 1 " and P in = Up 1 ™. For the vacuum 
input fields we define the standard quantum limit (SQL) 
to be: 

(Q^i-n)Q±m) = {pr(-^)p±m) = m (is) 

and all other correlations functions of the in fields are 
vanishing for a vacuum state input. Clearly A must de- 
pend on our choice of tp since for ip = Q± Lt describes the 
fluctuations of the individual signal and idler fields of the 
ND-OPO, whereas for ip ^ Q5_ nt describes inter-beam 
correlations. As we shall see below in the case of sym- 
metrical signal/idler losses and perfect inter-beam corre- 
lations we measure the strongest inter-beam correlations 
when choosing ip — tt/4. Hence we define A(ip = 0) = 1 
and X(ip ^ 0) = 2, since in the latter case we must 
compare our inter beam-correlations to a SQL given by 
the correlations between two independent coherent state 
fields. It is worth noting here that only very close to the 
ND-OPO threshold are the inter-beam quantum corre- 
lations perfect and the noise of the signal and idler can 
cancel each other exactly. However the single beam noise 
grows rapidly as the ND-OPO pump parameter e grows. 
Hence due to imperfect quantum correlations between 
signal and idler fields it turns out that since this growing 
noise cannot be canceled, it is more advantageous to mea- 
sure not the symmetric combination of the quadratures 
described by ip — 7r/4, but rather a somewhat scaled 
combination with ip = 7t/4 + <5(e)||. 

The correlations of the outgoing quadratures 

(0) and fjj) to be 



= A < I + Aen- 



2e 




( cos 2 ip \ 
V sin 2 ip J 


+W sin2 2 t) 

\ cos ip J 


+ (A 2 + I + e 2 )sin2^ 


(i) 




A 2 + (£ + A) 2 




A 2 + (E - A) 2 





(19) 



Here we have defined the normalized fluctuation fre- 
quency A = £!/^/7+7_, the signal/idler escape efficien- 
cies rj± — k±/j±, the generalized OPO escape efficiency 
for the signal and idler fields r\ = the escape 

efficiency asymmetry a = yJ~rj+~pqZ, the loss asymmetry 



function A = (p + p 1 )/2, with p = vT+Tt- an d the 
effective pump parameter E = \Je 2 + A 2 — 1. From ( fl9| ) 
we find for ip — the single beam noise spectra 

((^i 1 *) 2 ) = 1 + ? /± [A2 + (B+A)^f[AS + (B-A)S] ( 20 ) 
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FIG. 2: Correlation spectra for k+ 
asymmetry. 



and increasing loss 



FIG. 3: Correlation spectra with increasing 7+ for the opti- 
mized quadrature Q™*. Here tpo has been used. 



These are independent on 9, meaning that the individ- 
ual signal and idler fields contain phase insensitive ex- 
cess noise relative to the standard quantum limit here 
being 1. As the threshold is approached E — > A, and 
we find the noise to diverge at zero detuning. Naturally 
the linearized approach taken here breaks down in this 
limit, where we must take pump depletion into account 
as the ND-OPO gain becomes large and much energy 
is transferred from the pump to the fluctuations of the 
signal/idler fields. 

By setting ip = 7r/4 we find from ( |l9| ) the inter-beam 
correlations 



<(Q±T). /4 = 2{ 



1 ± Aerj 



A 2 + l + £ 2 ±((T+g~ 1 )£ 

[A2 + (B+A)2][A2 + (_E-A)2] 



(21) 

It is easy to check that in the case of symmetric losses 
p = a = 1, A=l and E = e, we recover from (^lj) the 
well known result for the cross correlations M 



<(Q° ± "*) 2 > 7r/4 = 2 H±n 



4c 



'A 2 + (c T l) 2 



(22) 



From ( p2| ) we find perfect inter-beam correlations at low 
frequencies as the threshold is approached e — ► 1, which 
reflects that tp — ir/4 corresponds to the natural sig- 
nal/idler linear combination to measure in the case of 
symmetric losses for the two modes. 

When allowing for non-equal losses of signal and idler 
modes, we see from Fig. |2j that a low frequency noise 
spike appears in the spectrum of fluctuations when 7+ is 
increased relatively to 7_. With p = 3 we find that the 
inter-beam correlations are virtually lost. This is because 
the simple sums and differences of signal idler fields, cor- 
responding to choosing tp — 7r/4, no longer are the best 
combinations to measure in order to observe strong inter- 
beam correlations. As we shall see by optimizing tp we 



can recover at least some of the inter beam correlations 
around a given frequency. 

A remark should be made about losses experienced 
by the signal and idler fields during propagation from 
the ND-OPO to the detectors. Any such losses includ- 
ing non-perfect detectors, limited local oscillator overlaps 
etc. for the signal/idler fields will reduce the quantum 
correlations on equal footing with non-unity escape effi- 
ciencies. Hence we define overall propagation efficiencies 
for the two field modes describing the survival probability 
for photons in these modes from birth in the process 
until detection. We denote these efficiencies £± for the 
two modes and for later reference we define them here as 



t- prop LO QE 



(23) 



Here we have defined the signal/idler propagation effi- 
ciencies r]j^ op , the overlap with the local oscillators of the 
heterodyne quadrature detections 77^° and the quantum 
efficiency of the detectors rf± E . In the realistic case of 
lossy propagation and non-perfect detection we should 
replace the signal/idler escape efficiencies in (|l9| ) with 
the overall efficiencies defined in (p3|). 

As mentioned above we can optimize the inter-beam 
correlations at a given set of parameters by choosing to 
measure the optimized linear combination described by 
ip — tpo, where ipo is found from 



tan 2ipo — 



1 



(24) 



Here Ao is the frequency at which we wish to minimize 
fluctuations. From ( p3J ) it is again clear that in the sym- 
metric limit where a = 1 there is no advantage of de- 
viating from measuring equally weighted gi"*(0+) and 
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FIG. 4: Correlation spectra now with fixed a and r\ and 
increasing 7+. This corresponds to the situation where k+ 
grows with the same rate as 7+. 



<7°"*(0_) observables corresponding to -00 = 7r /4- But in 
general we find a large variation of ipo with a for small 
values of A , hence there may be a big advantage of opti- 
mizing ip when signal/idler losses become non-symmetric. 
As Ao is increased the variation of ipo gets smaller since 
at higher frequencies the impact of measuring unevenly 
weighted signal/idler amplitudes gets negligible. 

By using ip — ipo hi (|l9|) and again varying 7+ relative 
to 7_ , we find the inter beam correlation spectra shown 
in Fig. |[ These have been optimized at Ao = 0. Here 
we see that the correlations at low frequencies persist 
despite the growing asymmetry of the losses, however 
since the escape efficiency 77+ decreases with increasing 
7+, we find the overall degree of inter-beam correlations 
to diminish when p deviates from 1. The importance 
of measuring scaled signal and idler fields is seen more 
clearly from Fig. [|, which corresponds to the somewhat 
unrealistic situation when p is increased, but a and r\ 
are left unchanged. Still keeping ip optimized we find no 
loss of inter beam correlations in an increasingly small 
frequency window around A = 0. Optimizing ip around 
a nonzero frequency is not possible without a reduction 
of correlations and in general the gain in correlations by 
optimizing ip is marginal when Ao > 1. 

As a conclusion to this section we remark that not only 
the quantum fields will show the behavior depicted in Fig. 
|l], but also classical fields. This means that if we inject a 
classical coherent field with frequency lu+ and phase </>+ 
into our ND-OPO, the nonlinear interaction will result in 
generation of an idler field of frequency cj_ and phase <p- 
such that (p+ + 4>- = Xj where \ is the pump field phase. 
This was utilized in the experiment to gain information 
about the pump phase of the ND-OPO. 



III. EXPERIMENTAL DETAILS 

We will now consider the experiment carried out to 
characterize the frequency ND-OPO described in the pre- 
vious section. This was divided into three parts: First we 
simply measured the squeezing in the degenerate mode of 
the OPO. Then we investigated the signal/idler cross cor- 
relations of the ND-OPO output and finally we recorded 
the signal/idler phases to confirm the relationships de- 
rived above and to visualize the EPR paradox for our 
kind of quantum mechanically entangled light. But be- 
fore considering the results we will just outline here some 
of the underlying experimental setup. 

This was made up of a single ended optical ring res- 
onator with a 10mm long Potassium Niobate crystal lo- 
cated in the smallest waist. Potassium Niobate is fa- 
mous for its high x nonlinearity which can be up to 
8 pm/volt. However it is equally notorious due to its 
nonlinear losses of which Blue Light Induced InfraRed 
Absorption (BLIIRA) is the most harmful (|, [|. In this 
experiment the nonlinear crystal had a single pass non- 
linearity of 0.015 ± 0.002 W _1 measured by second har- 
monic generation. The ND-OPO output coupler had a 
power transmission of 0.122 around the signal/idler wave- 
length of 860 nm. Residual losses in the ND-OPO were 
0.007 ± 0.001 excluding the BLIIRA. The cavity length 
was 81 cm resulting in a free spectrum range of 370 MHz. 
With a typical value of BLIIRA of 0.015 at 160 mW blue 
pump power we find the ND-OPO cavity bandwidth of 
8.5 MHz FWHM. 

The ND-OPO pump was produced by performing ex- 
ternal cavity second harmonic generation (SHG) of a 
Titanium-Sapphire (Ti:S) laser operating at 860nm. This 
was done in a single ended ring resonator also with a 10 
mm Potassium Niobate crystal placed in the waist. The 
cavity was 70 cm long and the input coupler transmitted 
0.078 at 860 nm. Residual internal losses were at the 
level of 0.009 again excluding BLIIRA. The cavity waist 
was 20 jUm. close to the optimum for nonlinear conversion 
of 16 pm |10|| . However since the resonator was pumped 
with up to 800 mW of light from the Ti:S a slightly larger 
waist was chosen to compensate thermal lensing in the 
nonlinear medium. From the SHG we got up to 300 mW 
of blue light (430 nm) for pumping the OPO. The mode 
matching of the SHG Gaussian mode to the ND-OPO 
mode was typically 0.95 measured by the reference res- 
onator shown in the setup on Fig. |^. 

The signal and idler fields emerging from the ND-OPO 
to be analyzed were separated in oscillation frequency by 
two free spectrum ranges, that is 740 MHz. In order 
to separate the fields spatially we sent them onto a two 
ended ring resonator with a length of 20 cm and coupler 
transmissions of 0.078 at 860 nm. The length of this fil- 
tering cavity was stabilized to maintain resonance with 
uj + , such that this frequency component was transmitted. 
The cavity linewidth was 38 MHz and the free spectrum 
range 1.5 GHz, meaning that less than 10 -3 of the idler 
field was transmitted when the filtering cavity was reso- 
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of RF spectrum analyzers. The quantum efficiencies of 
the photodiodes were measured to be higher than 0.98 
and the fringe visibilities on the 50/50 beam splitters 
were 0.97 and 0.98 for the lo- and oj+ fields respectively. 

As mentioned earlier we also investigated the degree 
of squeezing in the degenerate OPO output mode. For 
this purpose we inserted a mirror between the OPO and 
the filtering cavity. This mirror deflected the OPO out- 
put to a third heterodyne detector with a local oscillator 
extracted directly from the Ti:S laser. Hence the OPO 
output field was demodulated at exactly the frequency of 
degeneracy of the OPO. The fringe visibility in the de- 
generate heterodyne detector was 0.965 and the detector 
quantum efficiencies were also here higher than 0.98. The 
phase dependent quantum fluctuations were recorded by 
feeding the differential photocurrent of the detector pair 
directly into a RF spectrum analyzer. 



FIG. 5: The experimental setup. Long dashed lines are in- 
frared beams (ui± and lu), short dashed lines are blue beams 
(2lu) and full lines are electrical cables. Frequency shifted lo- 
cal oscillators for the lu± heterodyne detectors were generated 
with the acousto optical modulators AOM(±). For the role of 
the miscellaneous cavities see text. The final photocurrents 
were Fourier transformed and recorded in the RF spectrum 
analyzers SA(±). 



nant with ui + . Mode matching and internal losses of the 
filtering cavity limited the transmission of the signal field 
to 0.95. This cavity as well as the SHG and OPO cavities 
were locked using the Pound-Drever-Hall FM technique 
(□J. The mode cleaning cavities were locked to give a 
constant transmission at roughly 80% of their maximum 
transmission. 

From the filtering cavity the signal and idler fields 
propagated in different directions each to be heterodyned 
with a frequency degenerate strong local oscillator and 
detected on a pair of detectors. The local oscillators 
were produced by sending part of the Ti:S output in 
double pass through a pair of acousto optical modula- 
tors driven by a 185 MHz RF source. A small portion 
of the up-shifted light (a>+) was injected into the filter- 
ing cavity counter propagating to the ND-OPO idler field 
in order to stabilize the cavity length. Additional small 
portions of the ui± local oscillators were split off to have 
the possibility to inject them into the ND-OPO through 
a high reflector. The remainder of the local oscillators 
were passed through their respective mode cleaning cavi- 
ties to ensure a clean TEMoo mode and finally interfered 
with the ND-OPO outputs on separate 50/50 beam split- 
ters. The outputs from the beam splitters were detected 
on four Silicon pin diodes and the photocurrents from 
each pair of detectors we subtracted to make a balanced 
measurement of the quadrature phases of the signal/idler 
fields. The two differential photocurrents were led into 
a hybrid junction producing the sum and difference of 
these photocurrents to be Fourier transformed in a pair 



IV. RESULTS 

In order to learn about the ability of our OPO to pro- 
duce non classical quantum correlations we first analyzed 
the degenerate output mode. Here either increased or 
reduced quantum fluctuations should be observed de- 
pending on the relative phase difference between the 
heterodyne local oscillator and the OPO output. The 
phase of fluctuations in the latter was determined by the 
pump phase. Since the same mechanisms are responsi- 
ble for squeezing as well as signal/idler cross correlations 
the squeezing results yielded valuable information about 
what could be expected when observing the nondegen- 
erate OPO modes. Next employing the filtering cavity 
we analyzed the two nondegenerate modes closest to the 
degenerate mode. Here the cross correlations between 
signal and idler fields were measured while we tracked 
the local oscillator phases. This enabled us to visual- 
ize the EPR paradox now with frequency nondegenerate 
quantum correlated fields as opposed to the polarization 
nondegenerate fields of Ou et al [Q. We will discuss this 
in more detail below. 



A. Classical gain and squeezing 

When observing squeezing we pumped the OPO with 
160 mW of 430 nm light from the SHG stage. At this 
power we observed a BLIIRA level of O.Of 5±0.002 and by 
injecting a weak coherent field into the OPO via a high 
reflector we observed a maximum phase sensitive gain 
of 6.8. It can be shown |l^] that this gain is given by 
(1 — e)~ 2 , which gives us e = 0.62. Taking into account 
the blue mode matching of 0.95 this gives us an OPO 
threshold power of 395 mW. From ath putting 7+ = 7_ 
we find the power needed to drive the OPO above thresh- 
old to be 

{T + C + C P f 
Pth = 4E NL (25) 
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OPO escape efficiency. ?; 


0.847 




Heterodyning fringe visibility, ?7 i0 


(0.965) 2 
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Detector quantum efficiency, rf* E 


0.98 




Propagation efficiency, rf rov 


0.98 


6 


Overall detection efficiency, £ 


0.76 ±0.02 





TABLE I: The limitations to the degree of squeezing and 
their contributions. 



where T is the OPO coupler transmission, C are resid- 
ual OPO losses, Cp is BLIIRA and Enl is the single 
pass nonlinearity of the x medium. From our derived 
threshold power we now find E NL = 0.013 ± 0.001 W" 1 
in reasonable agreement with the value of 0.015 ± 0.002 
W _1 obtained from second harmonic generation in the 
OPO. 

Given the measured losses of the OPO, the quantum 
efficiencies of the detectors, the fringe visibility in the 
heterodyne detection scheme and measured propagation 
losses of 2%, we find using ( p3| ) an overall detection effi- 
ciency of quantum correlations of £ = 0.76 ± 0.02. The 
squeezing balance is summarized in Table @. The pho- 
tocurrent resulting from the balanced heterodyne mea- 
surement was sent into a RF spectrum analyzer to be 
Fourier analyzed in a 30 kHz frequency window centered 
at 1.15 MHz. Here the quantum noise of the OPO out- 
put (g 2 ) was recorded while the local oscillator phase was 
scanned at a rate of approximately 10 rad/sec. To es- 
tablish the standard quantum limit the OPO pump was 
blocked and the resulting vacuum input to the hetero- 
dyne detector was analyzed to give the vacuum noise 
level. By blocking all inputs to the heterodyne detec- 
tor the electronic noise floor was established. This was 
found to be 11.76 dB below the SQL when the local os- 
cillator power was 5 mW. A typical spectrum analyzer 
trace of the quantum noise is shown in Fig ^. From an- 
alyzing 10 such traces we found quantum fluctuations 
in the squeezed quadrature of the OPO output to be at 
— 4.0±0.1 dB relative to the SQL and in the anti squeezed 
quadrature 8.0 ± 0.2 dB of excess noise was found. Since 
our detectors had limited sensitivity the squeezing results 
were influenced by the electronic noise. After subtrac- 
tion of electronic noise we found the corrected squeezing 
of -4.5 ±0.1 dB and the anti squeezing of 8.3 ± 0.2 dB 
relative to the SQL. The degree of squeezing should be 
compared to the theoretical value found from |l3| 



4e 



A 2 ±(l±e) 2 



(26) 



a 

CO 

2 o 



a-2 
v 




0.1 0.2 0.3 0.4 
Time [sec] 



FIG. 6: The quantum fluctuations of the degenerate OPO 
mode. The trace was taken with the local oscillator phase 
scanned. A minimum of -4.0 dB squeezing relative to the 
SQL was observed. 



Loss source 


Signal 


Idler 


V± 


0.847 


0.847 


LO 
V± 


(0.98) 2 


(0.975) 2 


QE 

n± 

prop 


0.98 
0.91 


0.98 
0.96 


e± 


0.73 ± 0.02 


0.76 ±0.02 



Inserting our parameters in (|26|) we find an expected de- 
gree of squeezing of —5.2 ± 0.3 dB which we believe is in 
reasonable agreement with the observed value. 



TABLE II: The detection efficiencies for the ND-OPO signal 
(±) and idler (-) modes. 



B. Entangled state of two modes 

With our OPO characterized in terms of its ability to 
produce quadrature phase squeezed light we could now 
move on to investigate the quantum correlations between 
the frequency non degenerate OPO cavity modes adja- 
cent to the degenerate mode. Since these modes only 
were shifted ±370 MHz relative to the degenerate mode 
it was safe to assume equal output coupler transmis- 
sions and internal OPO losses for both modes. How- 
ever, since the signal and idler fields experienced dif- 
ferent optical paths, the detection efficiencies were not 
equal for the two modes, meaning that a ^ 1 in (|i~s|), 
whereas we had p = 1. The limiting factors on the in- 
ter beam correlations are listed in Table ||, from which 
we derive the parameters £ = '£+£_ = 0.74 ± 0.01 and 

a = V W£- = °- 98 ± °- 01 - B y usin S (H we nnd the 
optimum phase choice ipo = — 7r/4 + 0.0086 radians which 

inserted in ( fi9|) gives a predicted maximum degree of in- 
ter beam correlations of -4.9 dB. It should be mentioned 
that the correction to ip due to asymmetric losses was 
so small that we chose not to make any at all. This 
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FIG. 7: Quantum fluctuations of the signal/idler quadrature 
difference taken at ip = ipo and with various local oscillator 
phase settings, (a) Standard quantum limit for two indepen- 
dent coherent states. Everything is normalized to this level (2 
SQL), (b) Standard quantum level for a single coherent state. 

(c) The correlated observable (Q+ lt )i, g — q+(<p+) — q~(<p-). 

(d) Anti correlated observable (Q_ )<0 O ~ q+((j>+) + q-((j>-) ~ 
q+(<p+ + 7r/2) — q-(4>- + 7r/2). (e) <p kept constant while <f> 
was scanned (almost) linearly in time. The traces are taken 
at 1.15 MHz with a resolution bandwidth of 30 kHz. 

also applied to the correction of ip due to a finite de- 
gree of correlations. Hence within our uncertainties we 
used ip = 7r/4 throughout the experiment. The measured 
correlations are shown in Fig. |7], from which we find 
the noise reduction relative to two independent coherent 
state correlations of —3.8 ± 0.1 dB. Correcting again for 
electronic noise we get —4.3 ±0.3 dB noise reduction due 
to the inter beam quantum correlations. Again this is in 
reasonable agreement with our theoretical predictions. 
Fig. shows five traces: (a) The SQL obtained by block- 
ing the ND-OPO pump, hence only the vacuum input 
to the heterodyne detectors was being probed yielding 
the quantum fluctuations from two uncorrelated vacuum 
states, (b) The SQL from a single heterodyne detector. 
This level was found by blocking the local oscillator of the 
other heterodyne detector as well as the ND-OPO pump. 
Unblocking now all beams we obtained the trace (c). As 
documented by the trace, phase fluctuations in our setup 
were relatively small, so that the local oscillator phases 
could be maintained constant on a timescale of hundreds 
of milliseconds in a way that maximum inter-beam quan- 
tum correlations were observed. The same procedure 
was used to obtain the level of maximum anti correla- 
tions (d). Next one of the local oscillators phases was 
allowed to drift freely while the other phase was scanned 
at a high rate. This resulted in the oscillatory trace (e) 
showing the transition between observing the correlated 
<(Q+ Ut ) 2 ) and the anti correlated ((Q™*) 2 ) defined in 



&■ 

According to the criterion for state separability devel- 
oped in refs. |l4|, [ll| we must have 

( g+ (0 + )-M<M) 2 )+ (97] 

( g+ (0 + +7r/2) + (7 _(0_+7r/2)) 2 )<2 [Zi) 

in order to be able to claim state inseparability for the 
signal and idler fields and hence entanglement. Insert- 
ing the measured values of interbeam correlations we get 
0.37+ 0.37 = 0.74 < 2. Hence our signal and idler fields 
are indeed in an entangled quantum state. Since the 
right hand side of ( p^ ) sets the limit for separability we 
can quantify our degree of entanglement as -4.3 dB rel- 
ative to this limit. To our knowledge this is strongest 
entanglement demonstrated up to now using a narrow- 
band tunable lightsource. We are of course aware that 
no universal measure of the strength of entanglement has 
been established yet, hence our quantification above only 
applies for the Gaussian continuous variable case. 

In order to get a measure of the local oscillator phases 
we now injected a weak coherent field into the ND-OPO 
through a high reflector. The purpose of this exercise 
was to show that the sum of the local oscillator phases 
must be shifted by tt radians in order to switch from 
observing correlations to observing anti correlations as 
illustrated on Fig. [l]. The frequency of the injected field 
was w_, hence a coherent field with the frequency uj + 
was produced in the ND-OPO via difference frequency 
generation. Given that the injected field had the phase 
—9' relative to half the pump phase x/2 we know that 
the high frequency field was produced with the phase +9' 
relative to x/2- As a result we could by observing the 
DC interference fringe between these two fields and the 
respective local oscillators learn about the local oscillator 
phase settings relative to the ND-OPO pump phase. An 
example of the interference fringes is shown in Fig. |^. 
In order not to disturb the quantum noise measurements 
taking place in the heterodyne detectors the coherent am- 
plitudes of the injected field had to be small. Hence only 
about 100 nW hit the detectors and consequently the DC 
interference fringe visibility was only about 0.01 with a 
rather poor signal to noise ratio as a result. Again we 
left one local oscillator phase freely drifting but this time 
the other local oscillator phase was actively stabilized to 
produce either maximum or minimum inter beam quan- 
tum correlations. This was done using a standard dither 
and lock approach with an error signal derived from the 
spectrum analyzer output. Now the phase of the injected 
field was scanned linearly in time and the resulting inter- 
ferences with the local oscillators were recorded. Before 
commenting on the results displayed in Fig. |^ we should 
note that when 9' was varied in the negative direction for 
the ui- field, the change of 9' would be in the positive di- 
rection for the w+ field due to phase conservation in the 
difference frequency generation process. Hence we would 
expect the phase difference between the two interference 
fringes to be constant when the local oscillator phases 
were fixed to produce maximum or minimum quantum 
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FIG. 8: The DC interference fringes revealing the relationship 
between the local oscillator phases and the ND-OPO pump 
field. The left column corresponds to the situation when the 
local oscillator phases were locked to produce maximum anti 
correlations in (Q_)^ , while the right column are the results 
when maximum correlations were observed. The upper row 
are the raw data and the lower row are the local oscillator 
phases and their differences inferred from sinusoidal fits to 
the raw data. 



correlations between the two heterodyne detectors. If 
this phase difference was constant it would correspond 
to having the local oscillator phases adding up to a con- 
stant value as illustrated in Fig. M. We found the phase 
differences Og = 6± - 6'± and 9^ c = 0± - 6'± for the 
two frequencies in the cases of being locked to maximum 
correlations and maximum anti correlations respectively. 
This was done by fitting the somewhat noisy interference 
fringes to a sinusoidal function and the resulting phase 
variations as well as their differences are shown on the 
graphs of Fig. | We found 9% - 9 C = 0.53 ± 0.11 rad 
and 9 AC ~ 9 AC = -2.52 ± 0.06 rad, and we note that 
{6% - 9°) - {9 A c - 9 AC ) = 3.05 ± 0.13 rad. Ideally the 
last result should be ir since the sum of the local oscillator 
phases should be shifted by this amount in order observe 
the change from measuring a quiet to a noisy quadra- 
ture combination. As a final remark here we note that 
our setup was one big Mach-Zender interferometer from 
which the systematic phaseshift of 0.53 radians seen in 
99 — 9^_ was originating. The procedure described here 
for establishing the local oscillator phases was used in 
the following experiment, where the EPR paradox was 
revived for a short while. 



C. The EPR paradox revisited 

It is well known that the quantum correlated field 
quadratures of the ND-OPO output can be used to visu- 
alize the EPR paradox in its original form put forward 
by Einstein, Podolsky and Rosen in 1935 |H. This has 
been first demonstrated using a polarization non degen- 
erate OPO §. It is known that the non-negativity 
of the Wigner function describing the ND-OPO output 
prohibits us from violating local hidden variable theo- 
ries when measuring quadrature phase amplitudes [ Q |. 
However an experimental realization of the EPR paradox 
is still interesting especially seen in the light of recent 
proposals suggesting that Bell inequalities generalized to 
continuous variables can be violated by measuring parity 
like observables jl8[ |l^, . Moreover when considering 
quantum information transfer with continuous variables 
using the ND-OPO, like quantum teleportation || one 
must have a fidelity larger than 2/3 for outperforming 
classical systems consisting of linear optics and coherent 
states in terms of reproducing the input state pl| . This 
can only be achieved when more than 50% noise reduc- 
tion is present in the correlation between the signal and 
idler fields, which is just the condition required for real- 
izing the EPR paradox ||. We note that the 50% limit 
only applies when one uses iji — 7r/4. For other values of 
if) it turns out that the EPR paradox can be realized for a 
smaller degree of correlations. However demonstration of 
the EPR paradox with our choice of ip in our ND-OPO is 
equivalent to demonstration of its usefulness as a resource 
for continuous variable quantum information transfer like 
cryptography j2^, ^3|, teleportation j24j, storage etc. 

We refer to refs. |^6| and jl], || for the detailed dis- 
cussion of the EPR paradox. In brief we arrive at the 
paradox by inferring the idler quadratures <7_(0_) and 
+ 7r/2) from measuring the correlated signal vari- 
ables q + ((/>+) and q+(4>+ + tt/2). One can choose at will 
to measure either of the two variables, hence according 
to EPR they both should be "elements of objective real- 
ity" . However according to quantum mechanics the non 
commuting observables and <7_(0_ +7r/2) cannot 

have objective values simultaneously, hence the paradox, 
which is demonstrated if the product of the inferred un- 
certainties is ((g_ ^-)?) mf + 7r/2)) 2 ). n/ < 1. 

Turning now to the experiment we sent the photocur- 
rents from the two heterodyne detectors into a hybrid 
junction producing the sum and the difference of the pho- 
tocurrents in its two outputs. These outputs were fed into 
two RF spectrum analyzers in order to establish the noise 
density around 1.15 MHz. During this experiment one 
local oscillator phase was kept constant while the other 
was scanned linearly. Like in the previous section we also 
here injected a weak signal into the OPO in order to keep 
track of the local oscillator phases. Depending on the 
phase of scanned local oscillator the outputs of the hybrid 
junction yielded the observables q+(4>+ + 9 + ) ± g_ (</>_). 
That is, with 9 + = we measured the non commut- 
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<(q + (0 + +0 + )-q_(0j) 2 > 
12 




<(q + (0 +9)+q_(4>jy> 




1 2 
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FIG. 9: The noise densities in the quadrature difference (up- 
per trace) and sum (lower trace) obtained while the relative 
local oscillator phase was scanned. Both traces are normal- 
ized to the SQL of a single coherent state. Here the 3 dB 
lines indicate the SQL of two coherent states. 



ing observables q + (<fi + ) — q_ 

+ 1-{4>-) = Q+0+ 
the (+) junction (see Fig. 
measured — (<z+ (4>+) + q 
-( ff+ (0+)-g_(0_)) = -(g + 
in the (+) junction. From 
we measured a quiet observable in the 
and hence gained information about q_ (</>_ 




_) in the (-) junction and 
tt/2) - q_(0_ + tt/2) in 
When 9 + = 7T we 
in the (-) junction and 
+ + tt/2) +q- (</>-+ tt/2)) 
2|) we see that with 6 + = 
-) junction 
while with 



6 + = 7T we measured a quiet observable in the (+) 
junction and then learned about the conjugate idler ob- 
servable q~{4>- + tt/2). The result of a scan is shown 
in Fig. [| Here the fluctuations of the two measure- 
ments are shown normalized to the vacuum noise level 
of a single detector. After averaging over 12 pairs of 
traces like the ones displayed in Fig g we found the un- 
certainties ((q + ((/) + ) - q-((/)-)) 2 ) = -1.0 ±0.1 dB and 



((<?+(</>+ + tt/2) + g_f>_ + tt/2)) 2 ) = -0.9±0.1 dB both 
relative to the single beam SQL. These uncertainties re- 
sulted in an error of the same size in our estimation of 
the idler quadratures. Hence we could readily calculate 
the product of inferred variances in this experiment to be 
((q-(4>-)) 2 ) mf <(«-(*- + V2)) 2 ) m/ = 0.65 ± 0.02 < 1 
which was a clear demonstration of the paradox. 



V. SUMMARY 

We have demonstrated how the frequency non degener- 
ate optical parametric oscillator can be used to generate 
spatially separated entangled optical fields. In this paper 
we have used the OPO cavity resonances adjacent to the 
frequency degenerate mode, resulting in entangled fields 
separated by 740 MHz corresponding to two OPO cavity 
free spectrum ranges. In principle one could select any 
pair of modes placed symmetrically around the degener- 
ate mode, as long as the x' 2 ' phase matching condition is 
fulfilled for this pair of modes. The entanglement for our 
source is present over a bandwidth of approximately 14 
MHz FWHM depending on ND-OPO pump parameter 
e. This results in a high spectral brightness making the 
source well suited for addressing atomic transitions. 

In particular our source can be used for quantum mem- 
ory readout via interspecies teleportation of a collective 
atomic spin state onto the polarization state of beam 
of light |27|) . For this purpose two modes of light dis- 
playing entanglement in the Stokes parameters would be 
needed and only minor modifications to our source would 
be needed in order produce such states of light. 

We found a degree of squeezing of -4.0 dB of the fre- 
quency degenerate mode and we found inter beam cor- 
relations of the nondegenerate modes of -3.8 dB. Both 
results could be improved by approximately 0.5 dB if the 
thermal noise of our detectors could be suppressed. In 
this case the agreement with our theoretical predictions 
was reasonable. 

To gain information about the OPO pump phase we 
injected a weak coherent field with frequency w_ and 
thereby produced a oj + field via difference frequency gen- 
eration. These fields interfered with the local oscillators 
yielding us information about the local oscillator phases 
relative to the OPO pump field. Here we found that the 
sum of local oscillator phases changed by 3.05 ± 0.13 rad. 
when going from observing maximum inter beam corre- 
lation to observing maximum anti correlations. This was 
consistent with the theoretical value of tt. 

Finally we used the entangled fields of our ND- 
OPO to visualize the Einstein-Podolsky-Rosen para- 
dox in its original form with continuous quantum vari- 
ables. Here we found a product of inferred variances 
of ((g_(^)) 2 ) m/ <(^(^+7r/2)) 2 ) m/ = 0.65 ±0.02 
which shows the paradox since this is significantly smaller 
than the Heisenberg limit of unity. 
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